Surface and capillary transitions in an associating binary mixture model. 
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I. INTRODUCTION 

In the past few years there has been a considerable 
progress in the understanding of the physics of fluids 
in restricted geometries [J. The properties of Uquids 
or gases confined in narrow pores appear different from 
those in the bulk. A common trend of these systems is 
the fact that condensation occurs at pressures different 
from the saturation bulk pressure and the critical tem- 
perature for phase separation is generally lower than the 
corresponding bulk temperature. According to Nakan- 
ishi and Fisher y, U, |^ an incompressible fluid mixture 
(with isotropic interactions) under the effect of confine- 
ment displays a reduced tendency toward phase separa- 
tion at a fixed temperature, because the effective attrac- 
tive forces result weaker within a pore. The confinement 
also determines an effective dimensional reduction thus 
enhancing fiuctuation effects. In addition, the presence 
of adsorbing walls may induce a surface transition in a 
sub-critical fluid, even when this is in a single phase re- 
gion in the bulk. Such a behavior, which is now well 
understood by means of a mean field approach 0, 111 
is due essentially to the reduction of the contribution to 
the free energy from the the cohesive fiuid-fluid forces 
caused by the presence of the confining walls. There 
exists, however, a class of non-simple fluid mixtures for 
which the situation is different: these are the strongly 
associating mixtures. Strong directional attractive inter- 
actions such as hydrogen bonding or charge transfer com- 
plexing between molecules affect dramatically the prop- 
erties of the fluid or solid. An interesting phenomenon 
that occurs in such a mixture is the reentrant miscibility: 
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the mixture is completely miscible at low temperatures, 
up to a temperature where it separates into two liquid 
phases. As the temperature is further increased the com- 
plete miscibility reappears. The nicotine- water mixture 
is a paradigm for such a behavior jQj . The temperature- 
composition phase diagram shows closed-loops of immis- 
cibility, characterized by the existence of a lower and an 
upper critical solution point (LCSP and UCSP, respec- 
tively). Hirschfelder et al [lO| proposed a mechanism 
to explain the reentrant miscibility. The essential ingre- 
dient is the existence of highly directional strong attrac- 
tion (as the hydrogen-bond) between the unlike particles. 
The dispersive (isotropic) forces between these particles 
are assumed to be so weak as to favor segregation of the 
species. At low temperatures, the energy favors complete 
mixing. However, as the temperature increases, the en- 
ergy decrement due to the formation of directional bonds 
competes with the decrease of orientational entropy. As 
a consequence, the system phase separates when the tem- 
perature raises. This mechanism has been confirmed by 
theoretical studies on lattice 0, 0, 0| and continu- 
ous models , as well as in recent computer simulation 
studies [TsLllq. 

How the presence of substrates or the confinement in 
a pore affect the behavior of these mixtures (if in a man- 
ner similar to that of simple fluids, or if on the contrary 
these mixtures are peculiar) is an issue of great practi- 
cal and theoretical interest. The presence of selectively 
adsorbing substrates can interfere with the formation of 
bulk hydrogen-bonding, leading to a stabilization of the 
two-phase region with respect to the completely misci- 
ble phase. Computer simulation studies on confined lat- 
tice models ^3 show the appearance of surface transi- 
tions in a range of temperatures below the LCSP. The 
same mechanism has been invoked to explain some exper- 
iments 0, 01 that seemed to contradict the Nakanishi- 
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Fisher picture. 

In the present paper we shall study a lattice model 
of associating binar^y mixtures, introduced few years ago 
by Lin and Taylor _20, "211 and reconsidered recently by 
the present authors ,2^ 2^ 2Ac 2^ ■ Iii spite of its ide- 
alized nature, the model has the merit of allowing for 
an explicit analysis of its non trivial phase diagram. As 
an instance, the study of its bulk properties predicts the 
existence of closed-loops of immiscibility under constant 
pressure conditions. This analysis will be extended to 
the behavior of the mixture in presence of substrates and 
under confinement. 

The structure of the paper is organized as follows. In 
section II we briefly present the model and write explic- 
itly the equations for the order parameter profiles. Sec- 
tion III is devoted to analyze in detail the wetting behav- 
ior of the mixture in the presence of adsorbing substrates. 
In section IV we study the capillary behavior of the same 
system when confined between two parallel plates (a slit) 
and in section V we consider the same system confined 
in a network of very narrow cylinders. Capillary transi- 
tions (and surface transitions for the slab geometry) will 
be analyzed in terms of the physical parameters of the 
system. Finally, in section VI we present our conclusions. 
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FIG. 1: Schematic representation of the A-B LT model in 
the two-dimensional square lattice (the underlying lattice is 
not shown for clarity) . Relevant fluid-fluid and fluid-substrate 
interactions are also highlighted (see text for explanation). 



II. THE MODEL 

Some years ago, Lin and Taylor proposed a lattice 
model to describe a binary mixture in which the two un- 
like species can form highly directional bonds. In more 
detail, the Lin-Taylor (LT) model H is an A-B 
binary mixture lattice model, defined on a generic D- 
dimensional lattice of coordination number ly. The A 
particles are allowed to singly occupy any of the lattice 
cells. Every cell contains v equal sub-cells, one for each 
face of the mother cell. Again, each sub-cell can host 
at most one B particle, provided the mother cell is not 
occupied by an A particle. In addition to the repulsive 
interactions (which determine the previous occupation 
rules), one includes short-ranged interactions (see Fig. 

The interaction strength between a pair of nearest- 
neighbor A particles is eAA, <^ab between unlike species 
and Ebb between B-B particles. All these interactions 
are non- vanishing only if the particles share a face. 

The bulk properties, as well as the global phase dia- 
gra, rn, have been exhaustively analyzed in previous work 
[2^ I23L I23. I25I I. One of the main conclusions of such 
studies is that the appearance of closed loops of immis- 
cibility occurs when the interaction between A-A and 
A-B pairs are attractive {eAA < 0, e^s < 0), and 
eAB ^ i^AA + min(es_B, 0))/2. The latter condition is 
in complete agreement with the mechanism proposed by 
Hirschfelder et al (10| for the reentrant miscibility. In 
this work, we assume the bulk fluid-fluid parameters to 
have values: ^aa = —1 (it defines the energy scale), 
tAB = —0.65 and ebb = 0, so that we are in the regime 
where closed-loops of immiscibility appear. 



The presence of substrates introduces new interactions 
(see Fig. P). We shall assume that the substrates do 
not affect the fluid-fluid interactions, i.e. there is no sur- 
face enhancement. The interaction strength between an 
A particle and the substrate is Va, if the particle is in 
a cell in contact with the substrate, and zero otherwise. 
The interaction between the B particles and the substrate 
is orientationally dependent, i.e. its strength Vb is non 
vanishing only if the face of the B particle touches the 
substrate. The form of the substrate-fluid interactions 
are consistent with the bulk interactions. Hereafter, the 
substrate-fluid interactions will be assumed as attractive, 
i.e. Vk < and Vb < 0. For absolute values of Vb large 
enough, the substrate-B interactions compete with the 
bulk A-B interactions, so that one expect this type of 
substrate to enhance the demixing tendency for temper- 
atures lower than the corresponding LCSP. 



To study the LT model in presence of substrates, we 
consider the grand-canonical ensemble, where the system 
is in thermal and chemical equilibrium with a reservoir at 
flxed temperature T and chemical potentials ^a and ^b- 
The activities za and zb are deflned via za = exp(/3/ZA) 
and Zb = exp{Pi.iB), where (3 = 1/fcsT. After trac- 
ing over all conflgurations of the B particles, the system 
results to be isomorphic to a mono-component lattice 
gas, with renormalized interactions and chemical poten- 
tial. The details on this procedure can be found in Refs. 
[21I I2M I25I I . The grand-canonical partition function S of 
the LT model can be written in terms of an equivalent 
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Ising model as: 

S = {1 + 2zb + Zb 

( 1 + zbG-I^^^ 



exp (^N^H-"^ 



exp 



(1) 



where N is the total number of lattice cells, Ni, the num- 
ber of cells adjacent to the substrate, i'± the number of 
nearest-neighbors in the direction perpendicular to the 
substrate; is defined as the number of nearest neighbor 
cells in directions parallel to the substrate. Obviously, 
Finally, Zjj''^^{K,H,AHi) is the canon- 
ical partition function of the Ising model on the same 
lattice. Let us name K the effective coupling constant, 
H the effective bulk magnetic field, and AHi, the effec- 
tive surface magnetic field acting only on the boundary 
sites: 



^ ^ f3eAA 1 / y/l + 2zB + z%e-P- 
4 2^1 1 + zse-Z^'^s 



(2) 



and 715(1,3), the B particle density in the sub-cell s of 
the cell i is 



ns(i, s) 
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(8) 



if the sub-cell is by the substrate (see Fig. Otherwise, 
ns(z, s) is defined as: 
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where (siSj) is the 2-spin correlation function between 
the magnetization on the cell i and its nearest-neighbor 
j that it is associated to the sub-cell s. It is easy to see 
that in the homogeneous case these equations reduce to 
the expressions given in Refs. '2^, "231. 

The equivalence between the LT model and the Ising 
model, allows one to obtain the surface behavior of the 
mixture from the knowledge of the behavior of the latter. 
However, such a mapping, is not at all trivial and may 
lead in some cases to unexpected results. 



H = ilnz^-^^- J ln(f + 2^5 + ^16-/5--) 
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All equilibrium properties can be obtained from the 
knowledge of the partition function Eq. ^ by using 
standard thermodynamical relationships. In particular, 
the bulk pressure is obtained as: 
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lim fcsT = — — 
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where vq is the volume of lattice cell, and N = V/vq. 
The mole fraction Xyi(i) profile is obtained as: 



nA{i) 



(6) 



where n^(i), the A particle density in the cell i, is defined 
in terms of the local magnetization as: 
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III. WETTING TRANSITIONS IN THE LT 
MODEL 

Let us consider an A-B binary mixture in the presence 
of a planar, infinite and structureless substrate. Such a 
substrate confines the fluid to the half-space z > 0. When 
the fluid (e.g., the A-particle rich phase, a) is at bulk 
coexistence, the substrate can promote the formation of 
drops of the opposite coexisting phase (the B-particle 
rich phase P) on it, characterized by a contact angle 6. 
When 6 = 0, the drops spread over the substrate, i.e. 
the phase /3 wets completely the a-substrate interface. A 
finite value of 9, instead, corresponds to a partial wetting 
situation. The crossover between a partial to a complete 
wetting regime is called the wetting transition [2^. The 
surface tensions involved in the problem ((Tq.s, f/Ss and 
ctq/j, corresponding to the a-substrate, /9-substrate and 
a-f3 interfaces, respectively) are related to the contact 
angle via the Young equation [23| : 



O-as = Crps + CTap COS 9 



(10) 



Thermodynamically the surface tensions correspond to 
the excess grand-canonical free energy per interfacial 
unit area (respect to the bulk) [23|. Since both 
and aps have an analytical behavior at this thermody- 
namic state, the wetting transition is a (surface) ther- 
modynamic phase transition, as it can be seen from Eq. 
H10|l and the dependence of 9 on the temperature. As 
usual, appropiate order parameters corresponding to the 
wetting transition can be defined as first derivatives of 
the surface tension respect to thermodynamic fields. In 
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FIG. 2: Qualitative wetting phase diagrams for the LT model 
for Vb = 0. Thick solid lines correspond to critical lines, 
dashed lines to the A pure gas-liquid coexistence lines and 
thin solid lines to the wetting transition lines. The symbol 
s\a\l3 {s\l3\a) means the bulk conditions where there exist 
wetting of the /3-substrate (a-substrate) interface by the a 
{13) phase, respectively (see text). The bulk conditions under 
which there is partial wetting are also shown. 



simple fluids, the excess number of particles per inter- 
facial unit area or adsorption F (that corresponds to 
— (9(Tq.s/9^)t, being /j, the chemical potential) is the nat- 
ural order parameter choice. In mixtures, other choices 
are available, e.g. the volume integrated excess mole frac- 
tion per interfacial unit area. Microscopically, the wet- 
ting transition manifests itself as the growth of a macro- 
scopically thick [3 layer intruding in the a-substrate inter- 
face composition profile. Three different types of wetting 
transition may arise. If the layer width jumps suddenly 
from a microscopically finite value, the wetting transi- 
tion is first order, and an off-coexistence surface transi- 
tion (prewetting) is associated to it. If the layer width 
diverges continuously, the wetting transition is named 
critical. Finally, the wetting transition can occur as the 
accumulation point of an infinite sequence of first-order 
transitions in the layering wetting case (these transitions 
extend to the off-coexistence region). The conditions for 
such a transition depend subtly on the fluid-fluid and 
substrate-fluid interactions, and interfacial fluctuations 
can play an important role . 

In the LT model two-phase coexistence occurs when 
H = and K > Kc, where Kc is the lattice dependent 
Ising critical coupling. Moreover, the wetting behavior 
is controlled by the ratio AHi/iy±K , that is the ratio 
between the surface energy and the loss of bulk energy 
due to the presence of the surface. If such a param- 
eter is positive the substrate favors the a-phase to in- 
trude between the substrate and the /?-phase. On the 
other hand, if AHi/iy±K < the phase f3 intrudes be- 




FIG. 3: Wetting phase diagram for the square lattice LT 
model. The bulk parameters are fixed to eas = 0.65eAA < 0, 
eBB ~ 0, and Va = 0.25eAA- The thick solid line corre- 
sponds to the bulk critical line, and the thick dashed line to 
the liquid- vapor transition line of the A pure fluid. The thin 
solid lines correspond to the wetting transition lines for differ- 
ent values of Vb- Hereafter the reduced units are defined via 
T* = ksT/leAAl and p* = pvo/\eAA\, where vq is the volume 
of a lattice cell. See explanation in text. 



tween the substrate and the phase a. We shall focus on 
the coexistence side H = 0+ {H = for AHi < 
{AHi > 0), respectively (otherwise, only partial wetting 
can occur). If the absolute value of the parameter is 
larger than one {\AHi/i/±K\ > 1), then all coexistence 
surface corresponds to a complete wetting situation. If 
such a condition is not fulfilled, one will observe a par- 
tial wetting situation up to a value (dependent on 
lAHi/v^Kl) and beyond it a complete wetting situation 
up to = Kc- 

Let us consider first the case Vb — 0, when AHi de- 
pends only on T, regardless the value of zb- Moreover, 
the sign of AHi is uniquely determined by the differ- 
ence between Va and i'±eAA/'2 over the whole temper- 
ature range. One can envisage various cases (see Fig. 
I2J): if Va 0, the coexistence corresponds to a sit- 
uation of complete wetting of the interface substrate- 
a phase by phase /9, denoted by the symbol s|/?|q;. If 
i^±£aa/2 <Va<0, there exists a region of complete wet- 
ting s|/3|a and a lower pressure region of partial wetting. 
As Va decreases, the region of complete wetting s\(3\a 
shrinks and eventually disappears for Va — i'_\_eAA/'^, 
because AHi = 0. For Va < i'±£aa/'2, the sign of AHi 
changes and a region of complete wetting of the inter- 
face substrate-phase (3 by phase a appears (s|a|/?). As 
Va increases, the region of partial wetting reduces and 
disappears for Va > v±£aa- 

When Vb ^ 0, it is possible that in correspondence 
of different state points the substrate favors either the /3 
or the a phase. To study such an effect it is necessary 
to locate the coexistence states where AHi = 0. Since 
Vb and cab are negative, for a temperature T all the 
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FIG. 4: The same as Fig. Elbut for Va = 0.75eAA. The dotted FIG. 5: The same as Fig. glbut for Va = (-aa- See text for 
hnes correspond to the coexistence points at which A_ffi = 0. explanation. 
See text for explanation. 



coexistence states with pressure p larger than p„„, (T) = 
p(T, H = 0, AHi = 0) correspond to AHi > 0, whereas 
states with lower pressures correspond to AHi < 0. It is 
easy to prove from Eq. that a necessary condition for 
the existence of this line is that + < (Va/j^x) < 

caa/'^- However, AHi changes sign for a temperature 
T only if the value of the pressure which corresponds 
to AHi = is less than the critical pressure at such 
temperature. 

In order to illustrate such a phenomenology, we shall 
study the 2D square lattice case (v = A, v± = 1). In this 
case, only critical wetting can occur (fluctuation effects 
preclude any other possibility) and the value of at 
which the transition of critical wetting occurs is given by 
solution of the equation p9l | 

exp(2if„)[cosh(2i^„) -cosh(2AiJi)] = sinh(2is:„) (11) 

Let us first consider the case i^xeAA/2 < Va < 0: 
AHi < 0, for arbitrary values of T, zb y Vb- More- 
over, AHi becomes more negative as Vb decreases, at 
fixed T and zb (or pressure p). The typical diagrams 
of complete wetting are displayed in Fig. 13 The region 
between the wetting transition line and the critical line 
corresponds to a situation of wetting s|/3|q;, and the re- 
maining coexistence region corresponds to a situation of 
partial wetting. The complete wetting zone increases as 
Vb becomes more negative. 

When {Va/va.) < (aa/"^, the substrate also favors 
the adsorption of particles A. This determines a com- 
petition between A and B species. In general, when 
{Va/i^±) — ^aaI'^ < O.leAA /3-phase preferential adsorp- 
tion can only occur if Vb — {Va/v±) < cab — ^aa- Fig- 
ures ^ and El display the wetting diagrams for eAA < 
{Va/i^±) < eAA/'2 and {Va/v±) < eAA, respectively. In 
the first case, the complete wetting transition line orig- 
inates at T = and terminates at the temperature of 
complete wetting for the A-pure liquid in the interface 



substrate-A vapor. For values of Vb — 0, AHi > for 
all conditions of coexistence, and the configurations are 
s|a|/3. For Vb < (V4/i^±) + eAB — eAA there appear states 
characterized by AHi < 0. As a consequence, within the 
coexistence region at the left of the line p„m (T) the sub- 
strate adsorbs preferentially particles B, and the wetting 
is of type s|/3|q!. On the other hand, at the right of the 
curve Pnto(r) the wetting is of type s|a|/3. The line of 
wetting transition touches tangentially the critical curve 
at the point where the curve in turn pnw (T) crosses the 
critical line. Such a point moves monotonically toward 
higher temperatures as Vb decreases, driving with it the 
wetting transition line and the Pnw(T) curve. As also 
happened for the i^j_eAA/2 < Va < case, the wetting 
transition line from a s|/3|a complete situation to partial 
wetting moves to lower pressures, converging toward the 
A-pure liquid- vapor coexistence line for Vb —^—oo (how- 
ever, along this line there is complete wetting by liquid 
of the vapor-substrate interface over all the temperature 
range) . 

When Va < v^f-AA (see Fig. for values Vb near 
zero, the coexistence corresponds to complete wetting 
s|a|/3. Only when Vb < (Va/i'j.) -t- cab — caa states 
of partial wetting reappear, near the line Pnw{T). The 
behavior of the latter as well as the wetting transition line 
is similar to the one observed in the previous case, with 
the difference that the right branch of the wetting tran- 
sition line (that corresponds to a transition from partial 
to s|q;|/? complete wetting) begins at zero temperature. 

One expects that the behavior of the wetting transition 
for the three-dimensional case will be qualitatively simi- 
lar to the two-dimensional case. However, the complete 
wetting phase diagram for the three-dimensional Ising 
model is not totally understood yet. Theoretical 
and simulation studies [sj show, in the case of strictly 
short-ranged forces and in the absence of an enhanced 
surface spin-spin coupling, that the wetting transition is 
continuous only if is smaller than Kji, where is 



the value of the couphng constant corresponding to the 
roughening transition {Kji « l.S'oKc for the simple cu- 
bic lattice). However, if such condition is not fulfilled, 
the wetting transition will occur through a sequence of 
first order layering transitions. Theoretical studies have 
claimed that the wetting transition can be weakly first- 
order instead of second order [s^l- In the present work 
we will consider the Ising model in the framework of the 
mean-field Bragg- Williams approximation: 

lini -L In ZIs^ng {K, H, AHi) = AiJimi 

N±^co iV 



1 + rrij \ , / 1 + 
' In ' 



1 — m-j 



In 



1 — m. 



(12) 



where N± is the number of cells in a layer, an {rrij}, 
j — 1, 2, . . . oo is the equilibrium magnetization profile 
obtained by solving the Euler-Lagrange equations: 

nij = tanh(i^||/^mj + v^Klnij^i + rrij+i) + H) (13) 

where mg = AHi/v^K . This set of equations must 
be solved in conjunction with the asymptotic condition 
nij — > rrif, for j ^ oo, where rrib is the spontaneous mag- 
netization per site in the bulk case with coupling constant 
K and magnetic field H. In general, different solutions of 
Eq. (|13|l are found. The equilibrium profile is the solu- 
tion that maximizes the functional H12|) . This approach 
neglects capillary fluctuations, that it is equivalent to as- 
sume Kfi = Kc and the wetting transition occurs as a 
sequence of layering transitions |33l Is^ l . Once the equi- 
librium magnetization profile is computed, the surface 
tension a is obtained via: 



in + pV)/A 



(14) 



where = — /csTlnS, with S defined by the Eqs. ^ 
and (|12|l . and the pressure p by Eq. (O. The adsorption 
r, that it is the relevant order parameter in the wetting 
transition, is defined as: 



a) 



(15) 



where XaU) is the A molar fraction profile Eq. © by 
using the approximation (siSj) w mirrij. Finally, X\ is 
its bulk value. 

We have studied the simple cubic lattice, for which 
h' — 6 and vj^ ~ 1. The magnetization profiles have 
been obtained by an iterative method of the Eq. H13() 
for j = 1, . . . , Nmaxi imposing the condition rrij = mt 
for j > Nmax- In order to sample the solution space, 
we use sharp kink profiles {rrij = -t-1 (—1) for Ai/i > 




FIG. 6: Wetting phase diagram for the simple cubic (sc) LT 
model in the mean field approximation. The substrate-fluid 
interactions are Va = ^aa and Vb = 1.25eAA. The critical 
line (thick solid) , A-pure vapor-liquid transition line (dashed) , 
and p„w (T) line (dotted) are plotted. The wetting transition 
points (diamonds) are also represented. The line that joins 
them is only for eye-guide. 



(AiJi < 0), respectively, if j < jo, and rrij = rrih for 
j > jo) initial magnetization profiles. Starting the 
iteration with different initial profiles, in general, the al- 
gorithm converges toward different local minima of the 
free energy. As already mentioned, the global free energy 
minimum gives the true equilibrium magnetization pro- 
file. Finally, we have considered different values of Nmax, 
since the equilibrium profile in a complete wetting situ- 
ation has a strong size dependence. On the contrary, no 
appreciable finite width effects occur in a partial wet- 
ting situation for Nmax large enough. In order to locate 
the wetting transition, we studied the equilibrium profiles 
at bulk coexistence {H — 0) and constant temperature, 
varying the pressure. 

The wetting phase diagram for Va — ^aa and Vb ~ 
1.25eAA is plotted in Fig. As one can see, the quali- 
tative behavior of the wetting transition is similar to the 
one observed in the two dimensional case. However, the 
wetting transition is associated to series of first order lay- 
ering transitions between surface phases with the same 
surface tension, but different adsorptions. The behav- 
ior of the layering transitions at low temperatures shows 
surprising features. Fig. H shows the first four layering 
transitions for T* = ksT/leAAl = 0.25. It is observed 
that the first layering transition does occur for pressures 
higher than the bulk critical pressure. Consequently, for 
a fixed pressure between the first layering transition criti- 
cal point and the bulk critical point corresponding to this 
temperature, the layering transition will extend to lower 
temperatures than the critical temperature correspond- 
ing to the LCSP. This behavior is found only for s\(3\a 
complete wetting in the low temperature zone. The pres- 
sure is an increasing function on H for fixed temperature 
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FIG. 7: Layering transitions for T* — 0.25 for the sc LT 
model in the mean-field approach. The interaction couplings 
are the same as in Fig. El The dashed line correspond to the 
bulk coexistence, and the solid lines correspond to the bulk 
conditions corresponding to the layering transitions. In the 
inset, the coexistence adsorptions for the layering transition 
(only the first four transitions are plotted). The arrow shows 
the bulk critical pressure. 



T and zb (so both K and AHi are fixed). Consequently, 
as K > Kc for the layering transition critical point, the 
corresponding pressure can only be bigger than the bulk 
critical pressure if the transition occurs for H > 0. This 
implies that AHi must be negative, and thus the wetting 
is of type s|/3|q;. The physical interpretation of this phe- 
nomenon is that the surface adsorbs preferentially B par- 
ticles. Such a mechanism competes with the A-B bonding 
(at least close to the substrate). Therefore, the tendency 
to phase separate in a fluid layer close to the substrate 
is enhanced (because the A-B bonding favors mixing) for 
temperatures lower than the (lower) critical temperature. 
This fact was also observed in computer simulations of 
a different model of associating binary mixture in Ref. 
jl7j |. However, we must stress that this phenomenon is 
not related to an enhanced coupling (i.e. the interactions 
between particles in the first layer close to the substrate 
are increased respect the bulk values), that it is known 
to stabilize surface transitions to temperatures higher to 
the bulk critical one in the Ising model . On the other 
hand, the competition between substrate preferential ad- 
sorption and the A-B bonding should be also relevant for 
prewetting transitions. 



IV. CONFINEMENT BETWEEN SYMMETRIC 
WALLS 

We turn to consider the behavior of the Lin- Taylor 
model when confined between two identical parallel sub- 
strates. Confinement in the two-dimensional case corre- 
sponds to a quasi-one-dimensional situation. Thus fluc- 



tuations preclude the existence of capillary transitions 
if the interactions are short-ranged. Consequently, only 
the 3D case will be studied. We shall consider a sim- 
ple cubic lattice within the mean field approximation. 
The substrate-fluid interactions are set to Va — £aa and 
Vb ~ 1.25eAA, that stabilize the first (surface) layering 
transition at temperatures lower than the lower critical 
temperature. We also assume that the separation be- 
tween plates is No lattice units. The free energy of the 
confined mixture is given by Eq. , with the mean field 
partition function Zising (K, H, AHi ) given by 

lim In ZIs^ng (K, H, AHi ) = AHi {mi + tun^ ) 



No 

■E 

i=i 



-Km^ + iy±Kmjmj^i + Hr 



1 + rtij \ / 1 + m-, 
' In ' 



2 

1 — m., 



In 



2 

1 — m. 



(16) 



The equation for the equilibrium profile is given by 
Eq. (|13|l with the boundary condition N,,nax = Nu/2, 
for Nd even, or {No + l)/2, for No odd. The values 
of rrij for j > N^ax are obtained from the symmetry 
of the problem: i.e. nij = mNo-j+i- addition to 
the layering transitions (related to the phenomenology 
observed on the semi-infinite case), the capillary transi- 
tions results from the shift of the bulk coexistence under 
confinement. Near the center of the pore both phases 
have different magnetizations, in contrast with the layer- 
ing and prewetting transitions where the profiles differ in 
a localized surface region and they have the same values 
far enough from that region. In the Ising model, the shift 
is always toward higher values of K (i.e. lower effective 
temperature) and H negative {H positive) when AHi 
positive (A_ffi negative), respectively. 

The cases Nd = 1 and Njj = 2 can be studied ana- 
lytically since uij = m for all j. The bulk conditions for 
which the capillary transition takes place are shown in 
Fig. IHl In spite of the fact that the temperature range 
for which the capillary transition occurs is reduced, at 
low temperatures the pressure range is largely increased. 
Hence, at a given pressure, the LCSP shifts to lower tem- 
peratures and can even disappear. This means that for 
certain values of the pressure range there are no immis- 
cibility islands in the confined system, but the usual bell 
shaped coexistence line terminating in a UCSP. This is 
a direct consequence of the directional character of the 
substrate-particle B interaction. The kink that it is ob- 
served in the capillary critical line corresponds to the 
crossing to the Pnw{T) line. 

For larger values of Nd, the behavior of the capillary 
transitions are studied numerically. As N^ — > oo the 
capillary critical line must move toward the bulk critical 
curve. We shall consider the isotherm T* — 0.25, as a 




FIG. 8: Capillary transitions of the mean field sc LT model 
{Va = ^AA, Vb ~ 1.25eAA) for No ~ 1 (thin solid line) and 
Nn — 2 (dot-dashed line). In both cases, upper line is the 
capillary critical line, and lower line is the A-pure capillary 
gas-liquid transition. The bulk conditions under which there 
is a-/3 capillary coexistence are the regions enclosed by their 
respective lines, and the T = axis. By comparison, the 
bulk critical line (thick solid line) and the A pure gas-liquid 
transition line (thick dashed line). 
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FIG. 9: Xa — p phase diagrams for T* = 0.25 and different 
values of Nd: Nd = 1 (dashed-double dotted line), Nd = 2 
(dashed line), Nd = 5 (dotted line) and Nd = 15 (thick 
solid line), under the same conditions as in Fig. |H] The bulk 
transition (thin solid line) is also plotted for comparison. The 
inset corresponds to an enlargement of the zone around the 
bulk critical point. 



representative case of the low temperature regime. The 
excess free energy to reads: 

uj={VL+pV)/A (17) 

and converges to 2a when — > cxd, with a defined by 
Eq. (|14|l . The relevant order parameter for the capillary 
transitions is the average molar fraction XA^ defined as 

1 

i=i 

that converges to the bulk coexistence molar fractions as 
Nd oo. 

We have obtained the constant temperature phase di- 
agrams for T* = 0.25 (see Fig. ISJ. Various transitions 
are observed for No > 3. One of these transitions con- 
verges toward the bulk phase transition for large separa- 
tions and is identified with the true capillary transition. 
The remaining transitions are characterized by having 
similar average molar fractions and correspond to the 
layering transitions observed in the presence of a single 
substrate. The coverages (as defined by Eq. (|15|l ) cor- 
responding to the coexisting phases in these transitions 
converge very fast to the semi-infinite values, as expected 
from their surface nature. However, the number of these 
transitions depends strongly on Nd^ since the capillary 
transition, which depends sensibly on Nd^ competes with 
them [33, Hg- For instance, in the case of 3 < Nu < 10 
only the first layering transition appears, while the sec- 
ond appears for Nd > 10. 



The values of the capillary critical pressure for small 
Nd are larger than the corresponding bulk critical val- 
ues. By increasing Nd, the capillary critical pressure 
decreases up to a minimum, corresponding to a pressure 
less than the bulk critical pressure. From such value it 
converges monotonically to the bulk critical value (see 
figure I10|) . This fact means that along an isobaric the 
LCSP shifts to lower temperatures only for very narrow 
pores. Consequently, the stabilizing mechanism of the co- 
existence region below the bulk LCSP seems to be only 
effective for small values of Nd- The deviations of the 
capillary critical pressure with respect to the bulk value 
in the range 5 < Nd < 15 follow a power law Ap ^ N^^, 
with X = 1.04 ± 0.09. This behavior is consistent with 
the findings for the confined lattice gas for intermediate 
values of Nd reported in Ref. ■ The scaling behav- 
ior Ap ^ Njj^^" i^D^ i'^ the mean- field approach) pre- 
dicted by the Nakanishi-Fisher theory should be obeyed 
for larger values of Nd- However, our results do not 
extend to that region due to the numerical error in es- 
timating the capillary critical pressure close to the bulk 
critical point. Nevertheless, we do not expect a change 
in the tendency for large values of Nd and we expect our 
results to be qualitatively correct for all values of N p. I t 
is interesting to note that the slab approximation j35| . 
which corresponds to set rrij — m for all j, predicts that 
for sufficiently small T* the capillary critical pressure is 
larger than the bulk value for all the values of Nd- The 
failure of this approach is not surprising since it gives 
poor results for the confined lattice gas close to the cap- 
illary critical point as it overestimates both K and H at 
the capillary critical point [s^. This fact alters the bal- 
ance between the different mechanisms involved in the 




FIG. 10; Capillary reduced critical pressures p* versus No 
for T* = 0.25 (same conditions as Fig. |HJ. The dotted line 
corresponds to the bulk critical pressure. Inset: enlargement 
of the minimum zone. 



shift of the critical point, leading to the wrong result 
mentioned above. 
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FIG. 11: Schematic representation of the LT model adsorbed 
in a network of narrow infinitely long cylinders. The A parti- 
cles are represented as spherical particles, and the B particles 
as spherical particles with an interaction side (shadow circle). 
In this case, the coordination number is — 6. 



V. THE CYLINDRICAL PORE NETWORK 



volume reads: 



Finally, we consider the confinement of the LT mix- 
ture in a network of parallel infinitely long, cylindrical 
pores. This model can be understood as a crude approx- 
imation for adsorption in zeolites. It represents another 
instance in which the substrate can prevent the forma- 
tion of hydrogen bonds: the geometrical hindrance due 
to the confinement in very narrow nanopores. We shall 
assume that the pore diameter is so small that the A- 
B bonds can form only along the direction z parallel to 
the cylinder axis (see figure [TT|l . The cylinder diameter 
is taken to be less than twice the A-B collision diame- 
ter. In addition, the pores have diameters less than two 
A particle diameters and they form a two dimensional 
network with coordination number v±. We allow atten- 
uated interactions between particles in nearest neighbor 
pores, in order to allow a capillary transition. A similar 
zeolite model has been employed to study the capillary 
transition of methane adsorbed in ALPO4-5 37] . 

The in-pore fluid-fluid interactions have the same form 
as in bulk, with coupling constants eAA, ^ab and esB- 
However, the interactions between particles in nearest- 
neighbor pores are attenuated to eAA, ^ab and Ibb 
< although the A-B and B-B interactions 

directional character remains the same as before. This 
system behaves as a bulk lattice model with anisotropic 
interactions. The grand-canonical free energy per unit 



V 



2 



ln(l + 2zB + 4e-'^'««) + i/- 



(19) 



where Zb = ZBe~^ is the effective activity of the B 
particles oriented in the plane perpendicular to the cylin- 
der axis, and Zjsing is canonical partition function of an 
anisotropic Ising model characterized by couplings in 
the parallel directions to the cylinder axis and K± in the 
perpendicular directions, and an effective magnetic fleld 
H. has the same expression as the bulk coupling 
constant Eq. jSJ, while K± and H are defined as: 



^ -^ + lln Vi±li±5^ (20) 



H = H + vi 



\ £AA - f^AA 



2Va 



1 / l + 2zj3 + z|c-^^"^ 
4 ^ \ 1 + 2zb + z'^gC-P'-BB 



(21) 



with H given by Eq. (|3Jl. One can see that confinement 



shifts the effective magnetic field H with respect to its 
bulk value with a contribution analogous to that found 
for planar substrates (compare Eqs. (0)) and l(2U)- In 
addition, the confinement determines an anisotropy in 
the couplings. The A particle molar fraction Xa reads: 



HA 



(22) 



where ua = (l + m)/2uo is the A particle density, and 
(n^) is the density of B particles oriented in a perpendic- 
ular (parallel) direction to the cylinder axis, respectively: 



4wo 



zb 



+ 2Sb 
Zb + ' 



1 + ZBe-f^'^^ 

-P^AB 



zbB 



1 + 2Sb-\ 
dlnZistng/N 



dKi_ 



Zb 



i|e- 



■lie. 



1 + zbG 



-fie. 



2va 1 + 2Sb 



(23) 



ZB 



^2 f^ — pesB 



1 

2va V 1 + 2zB -r ^B 
1 / z_B + z|e " 



ZBe 



-/3eAL 



1 + ZBe-l^"^'^ 



zbg 



2vo \l + 2zB-\ 

d\nZising/N 

Zb t- ^b^ 



-B^ 



1 + ZBe-I^^^B 



vq1 + 2zb + zle-i^^BB ' 



(24) 



Phases characterized by different molar fractions Xa 
may coexist when the spontaneous (i.e. H — 0) magne- 
tization m = ±|m| 7^ 0. We consider the case eAA = — 1, 
eAB = —0.6 and ess = 0, that corresponds to a case that 
presents closed-loops of immiscibility in the bulk. On the 
other hand, as the hydrogen bond interactions are very 
short-ranged, we will take Iab = £bb = as physically 
sensible values. Moreover, eaa is set to be acAA, where 
a is a positive number less than one. 

In order to discriminate between different mechanisms, 
we chose the substrate-fluid interactions in such a way 
that H = H, i.e. Va = iy±{l - a)eAA/2 and Vb = 0. 
Hence, the preferential substrate adsorption effect (very 
similar to the one studied in the slab geometry) is sup- 
pressed and we can focus on the coupling anisotropy. 

We have studied the 2D square lattice case, for 
which analytical expressions for Zd ^, are known in the 
anisotropic case (38j . However, we expect that the fea- 
tures obtained for this case will be qualitatively correct 
also for the 3D case. From the exact solution, the val- 
ues of (i^_L, if||) that correspond to the phase coexistence 
satisfy the following relationship: 



sinh(2iC||)sinh(2ii:i) > 1 



(25) 



where equality holds only at the critical point. It is clear 
from such an expression that, even for (equal to the 
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FIG. 12: Capillary transitions in the network of cylindrical 
pores. The reservoir conditions corresponding to the capillary 
critical lines are represented for different values of a: a = 
0.2 (dashed line), a — 0.4 (dotted-dashed line) and a = 0.6 
(dotted line). The bulk critical line (solid line) and the A 
pure vapor-liquid line (long dashed line) are also plotted. In 
the inset, the phase diagrams corresponding to p* = 0.07 
for the same values of a (the symbol meaning is the same as 
previously) and their comparison with the bulk phase diagram 
(solid line). 



bulk value) less than the isotropic bulk value, the sys- 
tem can undergo a phase transition for K± large enough. 
Note that since K± = —f3aeAA/^, this condition is ful- 
filled in the low temperature region. This fact is also 
observed in the mean- field approach, for which the coex- 
istence condition is h'±K± + 2Kn > 1. 



Figure 1121 shows the reservoir conditions under which 
the capillary phase transition occurs for different val- 
ues of a. There also exist immiscibility islands for the 
capillary transition within the porous medium (see in- 
set in Fig. I12|l . but the corresponding LCSP occurs 
at temperatures lower than those of the bulk LCSP at 
the same pressure. As a, increases the temperature and 
pressure range for which one observes capillary transi- 
tion increases. Hence, the immiscibility islands become 
larger with a, but mainly in the high temperature region. 
In fact, the temperature corresponding to the LCSP is 
rather insensitive to the value of a, and it is always be- 
low the temperature corresponding to the bulk LCSP. 
Our findings can be understood by a simple physical ar- 
gument. The geometry forbids the formation of hydrogen 
bondings in the directions perpendicular to the cylinder 
axis. So, the confinement effectively reduces the A-B 
bonds respect to the bulk case, and consequently the 
phase separation is enhanced in the lower temperature 
regime, in complete agreement with our results. 
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VI. DISCUSSION AND CONCLUSIONS 

In this paper we have studied the effect on the phase 
equihbria of an associating binary mixture in the pres- 
ence of substrates and under confinement in very sim- 
ple geometries. We have focused on the conditions in 
which bulk system exhibits closed-loops of immiscibility. 
Our results show clearly that the presence of directional 
fluid-substrate interactions can stabilize purely surface 
transitions (such as layering and prewetting) and also for 
the phase separation under confinement in nanopores, for 
temperatures below the bulk LCSP. However, this mecha- 
nism becomes inefficient for the capillary transition when 
the surface-to- volume ratio is small, and consequently we 
doubt that it can explain the experimental results in Refs. 
0,^1, which would correspond to the latter case. The 
disagreement with the scaling predictions in the cases in 
which the temperature shift depends on the substrate 
remains, in any case, unsolved. The experimental data 
show a capillary critical temperature dependence on the 
film width that it is not consistent with the scaling pre- 
dictions ^ ^ and it could be a signature of an emerging 
new characteristic length scale. However, to our knowl- 
edgement there is no theoretical explanation of such a 
behavior. 



Some final remarks are pertinent. Firstly, the present 
model represents a simplified description of an associat- 
ing binary mixture, but keeps all the physical ingredients 
needed to reproduce the correct phenomenology. Conse- 
quently, we are confident that our conclusions can be 
relevant also for more sophisticated models. Secondly, 
fluid adsorption in nanotubes and zeolites is experimen- 
tally feasible and they could provide experimental probes 
of our theoretical predictions. 
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